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Abstract 

Recently several authors studied the segregation of particles for a system composed 
of mono-dispersed inelastic spheres contained in a box divided by a wall in the 
middle. The system exhibited a symmetry breaking leading to an overpopulation 
of particles in one side of the box. Here we study the segregation of a mixture 
of particles composed of inelastic hard spheres and fluidized by a vibrating wall. 
Our numerical simulations show a rich phenomenology: horizontal segregation and 
periodic behavior. We also propose an empirical system of ODE's representing the 
proportion of each type of particles and the segregation flux of particles. These 
equations reproduce the major features observed by the simulations. 
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1 Introduction 



Granular fluids are assemblies of moving macroscopic solid grains, includ- 
ing common materials such as; sand, powders, cereals, etc. Their ubiquity 
in various industrial processes (e.g., pharmaceutical, food) as well as in geo- 
physics and astrophysics motivates a better understanding of their properties 
(Q). Due to the inelastic interaction between the grains, these systems have 
a natural tendency to cool down by successive dissipative collisions leading 
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to a macroscopic state of rest. Nevertheless, when these systems are supplied 
with an external energy source (e.g., a vibrating wall), they reach a dynamic 
stationary state whose macroscopic behavior presents some similarities with 
a classical liquid or a gas. Moreover, granular fluids exhibit a rich and pe- 
culiar phenomenology for instance; non-gaussian velocity distributions, non- 
equipartition of energy, anomalous transport coefficients, a natural tendency 
to form dense and dilute regions (0), shearing instability, etc. Recently, a sig- 
nificant progress has been achieved in the understanding of these phenomena 
by applying the kinetic theory and the resulting hydrodynamical equations. 
Some years ago a generalized Maxwell demon experiment (3) has been devised 
that shed light on one of the most amazing properties of granular media. This 
experience consists of a box divided into two equal compartments by a ver- 
tical wall at the bottom and going up to the middle of the box, leaving the 
upper half of the box undivided. The box is filled with identical inelastic par- 
ticles submitted to gravitation and the energy is supplied by a vibrating wall 
(Fig. 1). For a high energy input, the system presents a homogeneous steady 
state, while when the energy input is decreased, a phase transition occurs 
and an asymmetric steady state prevails. This asymmetric steady state has 
been first studied theoretically by Eggers (0) in a slightly modified configu- 
ration with a separation wall extending throughout the entire height of the 
box. This wall contains a small hole at a given height through which parti- 
cles can be exchanged. In this formulation was assumed that the exchange of 
particles between the compartments corresponds to an effusive process and it 
was considered that the mean energy of the granular system decreases when 
the density increases. This last feature of granular systems is a consequence 
of the inelasticity of the collisions between particles and is responsible for the 
symmetry breaking in the Maxwell demon experiment. Indeed, an asymmet- 
ric fluctuation in the system induces a decrease of the mean energy in the 
denser compartment. Consequently, the escaping rate of grains from this com- 
partment is reduced and the grains have higher tendency to flow from the 
dilute compartment towards the dense compartment. This argument has been 
successfully verified bv both molecular dynamics simulations (5) and hydro- 
dynamic approaches (p). 

In this work, we focus on a Maxwell demon simulation with the same experi- 
mental setting as described in the first paragraph (Fig. 1) and characterized by 
two equal components. The binary system is composed by two kinds of smooth 
hard discs, whose dynamics depends on a large number of parameters. Here 
we have assumed that the particles are characterized by the diameters <Ti, a 2 
and by the same mass (m = 1). Concerning the restitution parameter a it was 
set equal to 0.95. The particles evolve under gravity, g (which acts along the 
y-direction) , in a rectangular two-dimensional box (Fig. 1) of height H and 
width W, and moving ballistically between successive collisions. The collision 
rule between grains i and j reads: 
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Fig. 1. (a),(b),(c),(d) Snapshots of the oscillation at different times. In (a), one 
plots the initial condition, with the small discs in box L and the large discs in box 
R. In (b), one may observe the transfer of small discs into box 2 and the vertical 
segregation. In (c), the transfer of large discs into box L has started. In (d), most 
of large discs are in box L. 

= Vi - ^(1 + a)e(e.vy) = Vj + ^(1 + a)e(e.v ii ) (1) 

with Vjj = Vj — Vj and = — r,-. The primed velocities are the post- 
collisional velocities and the unprimed ones the velocities before the collision, 
e = t^t is the unitary vector along the axis joining the centers of the two 
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colliding spheres. These rules have been often used to study fluidized granu- 
lar materials and have been proven to reproduce the specific phenomenology 
of granular fluids. The equations of motion of the N particles are solved nu- 
merically using the event-driven algorithm (0). In our simulations, the box is 
divided in two equal compartments, box L (left) and box R (right) of width 
W/2, by a wall parallel to the y-direction of height H/2. The walls are elastic 
except the ground which is a thermalizing wall at a given temperature Tq 
injecting the energy to the system. For this we used the dimensionless param- 
eter 7 = (-^) relating the gravity and the energy input. The system studied 
is composecf by 500 small particles and 100 big particles. The diameter ratios 
analyzed were 8.0 and 4.0 and in both cases we obtained the same dynamical 
scenario as we will show below. It is important to point out that this way of 
injecting the energy into the system has no intrinsic frequency and does not 
introduce any essential differences when compared to an injection of energy 
by a vibrating wall. 



2 Dynamical regimes 

In the one component case the dynamical system is described by the pro- 
portion of particles on the left n L (t) = — N ^ ■ It is known that when 7 is 
decreased a transition occurs (0). That is, for 7 > 7* on average the particles 
distribute themselves uniformly over the two compartments: n L (t) = 1/2 and 
for 7 < 7*, one of the compartments gets denser, leading to an asymmetric 
stationary solution n L (t) 7^ 1/2. This symmetry breaking is entirely due to 



the inelastic coupling between density and temperature explained in the intro- 
duction. For a binary granular fluid other key features come into play. Indeed, 
when granular media are shaken, particles of different sizes tend to separate 
vertically (@),(0). Moreover, the system exhibits non equipartition of energy 
(i.e., the mean energy of the different species is different (10)) and as we will 
see below this is reflected by a difference in the granular temperature for each 
kind of particle. The coupling of these phenomena with the horizontal spatial 
instability described above leads to a surprising new phenomenon, namely the 
emergence of an oscillatory regime. In the following, the state of the system 
is characterized by two dynamical variables: the proportions of small and big 
particles on the left (t) = 1 N and n^yt) = 2 N2 respectively. In our 
simulations we identified three well-defined regimes, separated by the critical 
values 7* and 7**, with 7* > 7** (Fig. 2): 

-Stationary homogeneous regime (7 > 7*^: If the shaking of particles is vigor- 
ous enough, the particles are uniformly distributed over the two compartments. 
In this regime, which has been also observed for a case of a monodispersed 
distribution of particles, the system homogenizes like a classical elastic fluid: 
nf(t)l/2 and n%(t)l/2. The stability of this solution has been corroborated 
by simulations with different initial conditions. 

-Periodic regime (7 < r y*)\ When 7 is decreased below the critical value 7*, i.e. 
when the energy input is lowered, a new behavior appears: the two densities of 
particles evolve periodically in time with the same frequency and separated by 
a phase delay (Fig. 3). This is what we call a granular clock. The mechanism 
behind these oscillations is the following. Suppose one starts the experiment 
with the small discs in box L and all the large discs in box R (Fig. l.a). 
Initially, the net flux of the small discs is directed towards the compartment 
of the large discs (Fig. l.b) since the latter suffer more collisions (the collision 
frequency is proportional to the diameter) and dissipate more energy, favouring 
the clustering instability. Progressively, the small discs accumulate in the box 
R and consequently a vertical segregation takes place. Our simulations show 
that the large discs are pushed upwards in the box (Fig. l.c). This vertical 
asymmetry, in the box R, increases the probability of escaping of the large 
discs and diminishes the probability of escaping of the small discs. This will 
induce a flow of large discs towards the box L, while small discs continuing 
to flow from L towards R. Hence, most of the large discs finish in box L and 
box R becomes rich on small discs (Fig. l.d). This situation is the opposite of 
the initial condition. Then, the whole process reverses and the initial state is 
recovered by the same mechanism. Hence, one oscillation has been taken place. 
Obviously, the motion of large discs drives the oscillation and this may explain 
the observed phase delay between the two signals. The power spectrum for long 
time runs of this regime consists of a well-defined peak (Fig. 2) surrounded 
by smaller peaks. When 7 is decreased below the threshold 7*, the size of the 
small peaks increases with respect to the main peak. The fact that this effect 
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Fig. 2. In the left figure, bifurcation diagram obtained from event-driven simula- 
tions, as a function of the energy parameter 7. The amplitudes are the oscillation 
amplitude, characterizing the oscillation regime and an asymmetry amplitude, de- 
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right figure, discrete Fourier Transform of the autocorrelation of the signal in the 
periodic regime (7 = 0.48) calculated for long runs (5 x 10 8 collisions). 

increases with the distance to the bifurcation point seems to indicate that 
some non-resonant intrinsic nonlinearities start to play a relevant role. This 
property has been also observed in the solutions of the model proposed in the 
next section. Finally, it is important to point out that in our simulations the 
mean energy of the large discs is lower than the small ones, as reported by 
Wildman and Parker for a binary system of grains (fiolb 

-Stationary asymmetric segregated regime (j < 7**j: Finally, when the input 
energy is lowered further a stationary asymmetric state appears characterized 
by a horizontal symmetry breaking. Initially, the system behaves as in the 
periodic regime but in this case the input energy is reduced and the number 
of large discs which escape from the box R is not large enough to launch 
the oscillation. As long as the vertical segregation takes place, small discs get 
trapped below the large discs in box R, in such a way that only large discs 
can flow from box R to box L. Hence, an asymptotic steady state is reached 
in wich box L is mainly composed of large discs. It is interesting to mention 
that Barrat and Trizac (II II ) have already observed such stationary asymmetric 
state for a binary Maxwell demon simulation without gravity. 



The succession of the above regimes is summarized in the bifurcation dia- 
gram of Fig. 2. In addition to the regimes described above we also observed 
a transition regime characterized by intermittent oscillations between the two 
opposite steady asymmetric steady states. This regime, which appears near 
the transition between the periodic and the asymmetric states, may be asso- 
ciated to the action of the internal noise on the dynamics of the macroscopic 
variables near the phase transition. One should note that the intermittent 
regime complicates the study of the second transition and does not allow a 
clear determination of the system frequency at that point. Nonetheless, re- 
sults obtained from the theoretical model proposed in the next section seem 
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Fig. 3. Typical time evolution of the densities in the periodic regime respectively 
for the simulations (7 = 0.58) and for the proposed model (A=0.1, R=3.0, A=1.67). 

to indicate that the frequency at the transition point goes abruptly to zero. 



3 Model 

In this section we introduce a simple empirical system of dynamical equation 
that reproduces qualitatively well the regimes observed in the above described 
simulations. It formulation has been inspired by a model proposed by J.Eggers 
(0) for the Maxwell Demon experiment with only one of disc diameter com- 
position. By using a hydrodynamic approach Eggers was able to find, under 
certain approximations, an expression for the escaping rate from one compart- 
ment in terms of the proportion of particles in the considered compartment 
only. If n denotes the proportion of particles, in say the left compartment, then 
its time evolution can be written as: |^ = —nF(n) + (1 — n)F(l — n) , where 
F(n), the escape rate per particle, is a decreasing function for large values 
of n. This feature represents the typical granular effect that in high density 
regions, the average kinetic energy is reduced by the high number of inelastic 
collisions, leading to a reduction of the escaping rate. In the right hand side 
appears a loss term, the escaping flux, and a gain term which has the same 
functional structure as the gain term but in this case it depends on (1-n), the 
proportion of particles in the right compartment. This last point guarantees 
the existence of a steady homogeneous solution: n — 1/2. In fact, the steady 
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state condition leads to three stationary solutions: the homogeneous one and 
two asymmetric configurations that are in good quantitative agreement with 
numerical simulations and these results are robust for a large class of functions 
F(n). 

Concerning the two-component problem described in the last section, we are 
looking for a model, as simple as possible, that can be able to reproduce the 
sequence of regimes obtained when decreasing the energy input, namely the 
homogeneous steady state, the oscillating regime and the asymmetric steady 
state. Given the inherent complexity of the system, we will only require in this 
work a qualitative agreement with the simulations. For doing this we define 
two variables in the left compartment; the proportions of small, n, and large, 
m particles. The corresponding proportion in the right half-box are (1 — n) and 
(1 — m) respectively. Here we assumed that the escape rate of the grains from 

the left half-box is of the form, F(n,m) = e~~ ^ Where A is a positive 
parameter representing the balance between energy input and inelasticity of 
collisions. This parameter decreases with the inelasticity (A — oo, for elastic 
collisions). The positive parameter R describes the fact that the cross section 
of the big particles is larger than for the small ones ( the cooling effect of the 
big particles is larger). 

The vertical segregation phenomena is taken into account, qualitatively, by 
introducing two more dependent variables and sr. They describe the in- 
fluence of segregation on the escaping rate in the left and right compartments 
respectively. We assume that these variables change the escape rate per parti- 
cle say in box L, by reducing this rate for the small particles, and by increasing 
it for the large particles, in agreement with the discussion of the previous sec- 
tion. By construction, the segregation variables should vanish when there is no 
vertical segregation. The next step consists in assuming that these variables 
reach respectively the values nm and (1 — n)(l — m) asymptotically. This ex- 
pression qualitatively reproduces the fact that there should be no segregation 
if a box is composed by only one component (n = or n = 0), while segrega- 
tion effects should be maximum if there are many small discs and large discs 
in that box, for instance (n = 1 and m — 1). Before going further, one should 
insist on the fact that the two supplementary variables Sl and sr have no di- 
rect physical interpretation, they are formal variables used in order to mimic 
the effects of segregation on the system dynamics. However, these variables are 
necessary in order to reproduce the observed phenomenology. Indeed, a model 
where sl and sr are not included does not lead to an oscillating regime. By 
considering all the points described in this paragraph the modeling equations 
can be written as: 



= - n F(n, m)(l - s L ) + (1 - n)F(l - n, 1 - m)(l - s R ) 
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Fig. 4. (a) Bifurcation diagram for the modeling equations (eq. [4]) , for A = 0.1, 
R=3.0. (b) Phase diagram, A = 0.1. 
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— = -mF(n, m)(l + s L ) + (1 - m)F{\ - n, 1 - m)(l + s R ) 
L = — A(nm — Si) — - = — A((l — n)(l — m) — sr) 



dt 



dt 



(2) 



The parameter A in the two last equation introduces the relaxation time of the 
segregation mechanism. By construction, the symmetric regime (n=|, m=|) 
is always a solution of the dynamics of the system. The dynamical equations 
2 have been solved numerically and they exhibit the three regimes observed 
by our even driven simulation. In figure 4 is presented the bifurcation diagram 
obtained for the model which as can be seen is very similar to the one obtained 
by even driven simulations (fig. 3). It is important to point out that other 
functional dependences for the escaping rate of the discs have been tested 
and indeed they did not show qualitative differences. For complementarity, 
in our analysis we introduced an additive noise into equations 2 and some 
features of the simulations that were not observed in the deterministic model 
are accessed. Namely, the intermittent regime and the presence of satellite 
peaks around the frequency of the deterministic oscillations. One must keep 
in mind that, although the proposed model is relatively simple, it has been 
formulated by adapting the escaping rate proposed by J. Eggers for a binary 
mixture of disc. 



4 Conclusions 



A Maxwell Demon experiment composed by a binary mixture of grains with 
two different diameters has been studied by means of even-driven simulations. 
It has been found that the phenomenology of this system is richer than the 
case of single diameter grain case. In the former case we found a symmetric 
steady state, an oscillatory regime, an intermittent behavior and a stationary 
asymmetric regime featuring grain segregation. We have performed simula- 
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tions over a wide range of parameters values for constructing the bifurcation 
diagram and for deducing accurately the power spectrum of the system. The 
observed oscillation phenomenon is mainly due to the combination of two 
typical granular effects: the vertical segregation with large particles above the 
smaller ones and the horizontal instability related to the clustering effect. 
Moreover, we have proposed a qualitative model for studying binary mixtures 
of grain inspired on the model proposed in (0) for the experiment with only 
one species of discs. In order to recover the phenomenology obtained by the 
simulations, we have introduced into our proposition a segregation mechanism 
through two new variables. The numerical integration and the linear stabil- 
ity analysis of our dynamical model (Eq. 2) showed its excellent ability to 
reproduce the major dynamical features obtained by simulations, namely; a) 
transition from a symmetric regime to an oscillation regime (Hopf bifurcation) 
characterized by a delay between the signal of the small and the large discs, 
b) passage from the oscillation regime to an asymmetric steady regime. 

It worthwhile to recall that binary demon experiments have been already per- 
formed by Van der Weele et al. (p). Even when these authors do not observed 
an oscillatory segregation regime (granular clock) they suggested, based on 
their experimental results, that the key factor point leading to an horizontal 
segregation instability may be the difference on the disk radii. In our simula- 
tions the radius ratio used was relatively large (8.0 and 4.0). This last may 
lead to think that the observed dynamics may be restricted to a particular set 
of parameters. Nonetheless, both our simulations and our modeling equations 
reproduce undoubtedly the major dynamical features observed by experiments 
and the results reported by A. Barrat and E. Trizac for a binary mixture of 
grai ns (without gravity). Recent numerical simulations developed by Marconi 
(|l3l ) seem to corroborate the existence of the granular clock existence. Finally, 
even when the proposed empirical equations (Eq.2) are relatively simple they 
undoubtedly reproduce the major intriguing features observed experimentally 
and by simulations. Presently, we are working on a re-formulation of our sys- 
tem of equations (Eq.2) which will allow to represent both the horizontal and 
vertical segregation in such a way that it permits to measure experimentally 
the parameters associated to the segregation mechanisms. The results of the 
last approach will be presented in a forthcoming paper. 
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